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Abstract 

A reflection model with three components, a specular spike, a specular lobe 
and a diffuse lobe is discussed. This model was successfully applied to de- 
scribe reflection of xenon scintillation light (A = 175 nm) by PTFE and other 
fluoropolymers and can be used for Monte Carlo simulation and analysis of 
scintillation detectors. The measured data favors a Trowbridge-Reitz distri- 
bution function of ellipsoidal micro- surfaces. The intensity of the coherent 
reflection increases with increasing angle of incidence, as expected, since the 
surface appears smoother at grazing angles. The total reflectance obtained 
for PTFE is about 70% for VUV light at normal incidence in vacuum and 
estimated to be up to 100% in contact with liquid xenon. 
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1. Introduction 

The scintillation light emitted in vacuum ultra-violet (VUV) region by 
noble gases/liquids (namely xenon, neon and argon) is used in particle physics 
detectors to detect very small energy depositions It is therefore desirable 
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to know in detail the VUV reflectance inside these systems. Moreover, this 
knowledge is also relevant to applications of VUV sources in synchrotron 
ultraviolet radiation facilities. 

Owing to a very low outgassing, chemical stability and high reflectance, 
PTFE (polytetrafluoroethylene) is a common choice for the inner walls of 
VUV systems j2[,@]. However, PTFE is an inhomogeneous dielectric whose 
reflection pattern is not sufficiently known in the ultraviolet. In this work we 
study the reflectance distributions of some fTuoropolymers, namely PTFE, 
ETFE, FEP and PEA. Most of the existing reflection distribution studies 
have been developed in the framework of computer graphics and computer 
vision and tested only for visible light [ij]. Experimental reflection data in the 
VUV region is difficult to obtain due to the light absorption by the oxygen 
molecules in air. We measured the reflectance of these fluoropolymers with 
xenon scintillation light (175 nm), using a specially built scatterometer, in 
which the atmospheric air can be replaced by a non absorbing gas such as pure 
nitrogen or argon j^]. Here we show data on the reflection by these surfaces 
along with the results of a model that describes effectively the observations. 

The roughness of the surface, which has effect on both diffuse and specular 
reflections, is described as an ensemble of micro-surfaces randomly oriented 
in space, followinga certain probability distribution function related to the 
surface structure |6(. A geometrical attenuation factor accounts for shadow- 
ing and masking by different micro- surfaces, especially at low grazing angles. 

We conclude that three parameters suffice to describe the observed re- 
flectance of unpolished samples of the aforementioned fluoropolymers, namely: 
i) the index of refraction, ii) the albedo of the surface and iii) a roughness 
parameter [3]. This model has been validated with data of various wave- 
lengths, including VUV light from the scintillation of xenon (175 nm) and 
light from UV and green LEDs. The reflectance distributions of the samples 
were measured both in and out of the plane of incidence. It was found that 
smoother surfaces are better described if a specular reflection spike is also 
included to account for coherent reflection at the surface average plane. 

Here we present a reflectance model of rough surfaces, considering that 
the light reflection is both diffuse and specular, with three components: i) 
the albedo due to internal sub-surface scattering of light, ii) the reflection at 
the surface by a myriad of random faces and iii) the coherent reflection spike. 

We conclude that although the reflectance of the PTFE is mainly diffuse, 
the intensity of both the specular lobe and the coherent spike can be signif- 
icant in the VUV, specially at grazing angles. Hence, this effect should be 
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taken into account in detailed analyses of the light collection in scintillation 
detectors and other systems. 

This work aimed primarily at describing the reflectance distributions of 
dielectric surfaces, namely in the VUV region, in a way suitable to be used 
in simulations and analysis of the light propagation through the volume of 
a scintillation chamber, up to the point of detection or absorption. In view 
of these applications this model was added as a class of the Geant4 Monte 
Carlo simulation toolkit [s| with the effect of improving the quality of the 
the simulations of experiments that rely on detection of the VUV light. 



2. Modelling the Reflectance 

The reflection pattern produced by a surface can be studied using the 
bidirectional reflectance intensity distribution function (BRIDF). This func- 
tion is defined as the ratio between the intensity reflected along a viewing 
direction, v, and the radiation flux, $j, incident at the surface along the 
direction i [9|] 

<"*M |1/srl (1) 

where d$ r /df2 is the flux reflected per solid angle towards the viewing direc- 
tion (see figured]). 

The reflectance distribution was modeled by a function g comprising three 
different contributions: the diffuse lobe go, the specular lobe g$ and the 
coherent specular spike g c , 

q = gv + Qs + Qc (2) 

In a dielectric the diffuse lobe is associated to internal scattering of the light 
that penetrates into the material and is scattered by sub-surface inhomo- 
geneities back into the incoming medium. For smooth surfaces this process 
is approximately lambertian, corrected by the Fresnel equations of reflection 

eD = ^caa9 r {l-F(6 u n/n )} {1 - F (0 r , n /n)} (3) 

7T 

where F is Fresnel reflection coefficient , no and n are the indices of refraction 
above an below the surface and pi is the albedo of the surface, i.e. the fraction 
of the light which is not absorbed during the internal scattering. 

The surface roughness, though small, is important to the specular re- 
flected light, including the coherently reflected spike. To describe the lobe of 
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Figure 1: The system of coordinates and directions. The vector i represents the direction 
of incident photons, v is the viewing direction. Primed angles are measured relatively to 
the normal to the local micro-surface, n', at the scattering point. 



specular reflection we adopt a geometric reflection approach in which spec- 
ular reflection occurs at a myriad of micro-surfaces randomly oriented in 
space, with some probability distribution function, P |6[ . The corresponding 
BRIDF function can be put in the form 

(l-C)FPG 
Ss = 4 cos 9 , (4) 

where G is a geometrical factor to account for the effects of shadowing be- 
tween micro-surfaces and C is a weighting factor corresponding to the frac- 
tion of light undergoing coherent reflection. The factor 1/4 cos 6^ relates the 
geometry of the solid angle made by the normals of micro-surfaces with the 
solid angle to where they reflect [z] 11]. The above process is in competition 



with the coherent reflection responsible for a specular spike in the reflectance 
distribution of the surface. This is a pure wave phenomenon and as such 



cannot be predicted by a geometrical approach to the reflection [12j]. The 
coherent contribution can be written as 

Qc = CFG (5) 
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where C measures the importance of the coherent reflection. 

The coherent reflection wave can be predicted by the model of Beckmann- 
Spizzichino. However, this calculation uses conductor boundary conditions at 
the surface, which may not be satisfactory on a dielectric interface. Notwith- 
standing, this model predicts that the intensity of the specular spike behave 
as an exponential function of the surface roughness and of the incident and 
reflected angles [13], as 



C = exp 



47T 

— ah (cos Oi + cos 9 r ) 
A 



(6) 



where 9i and 9 r are the incident and reflection angles respectively, cr^ is 
the standard deviation of the distribution of the height irregularities of the 
surface, assumed Gaussian, and A is the wavelength of the radiation. Thus, at 
large angles of incidence the surface appears smoother and the importance 
of the specular spike increases. However, the fits to the reflectance of the 
measured fluoropolymer dielectrics clearly suggest the following dependency 

f K 1 

C = exp < — — (cos 6^ + cos9 r ) > (7) 

where K is a constant to be fitted to the data. 

As stated above, the shape of the specular lobe is related with the surface 
structure function. The surface is modeled by an ensemble of micro-surfaces 
oriented at random whose normals make angles a with the macroscopic nor- 
mal and are distributed according to some (usually unknown) probability dis- 
tribution function P(a) (see figured]). We compared the predictions based on 
two possible probability distribution functions: i) the Cook- Torrance func- 
tion for V shaped micro- surfaces fig . Per, which is approximately gaussian 
for small roughness, and ii) the function deduced by Trowbridge- Reitz for 



ellipsoidal micro-surfaces [141 ]. Ptr- 



PT R (a-n) = - (8) 

7r cos a + tarn a) 

P C T(a;m) = „ 1 . exp (- tan 2 a/m) (9) 

i\m A cos 4 a 

The two functions are normalized to unity and the parameters 7 and m 
measure, in each case, the roughness of the surface. The fits to the measured 
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Sample 




n 


PL 


7 


K 


Extruded (_ 


L) PTFE 


1.50±0.03 


0.69±0.07 


0.055±0.007 


3.0±0.3 


Extruded ( 


) PTFE 


1.46±0.04 


0.63±0.07 


0.066±0.008 


4.3±0.5 


Pressed PTFE 


1.45±0.04 


0.74±0.07 


0.049±0.015 


1.7±0.2 


PFA 




1.44±0.04 


0.69±0.05 


0.057±0.006 


2.4±0.4 


FEP 




1.41±0.02 


0.22±0.04 


0.052±0.009 


1.2±0.4 


ETFE 




1.44±0.03 


0.13±0.01 


0.040±0.010 


1.0±0.2 



Table 1: Fitted values of n, p^, 7 and K for the samples measured at 175 nm. All these 
surfaces have been polished prior to measurement. The "Extruded _L" and "Extruded ||" 
refer to cuts perpendicular and parallel to the extrusion direction. 



data suggest that the Trowbridge- Reitz function is more suitable to represent 
the probability distribution of the micro-surfaces (see figure Ej). Therefore, 
in what follows we adopt this probability function unless otherwise stated. 

The geometrical factor G mentioned above, that corrects for shadowing 
and masking effects by adjacent micro-facets, was borrowed from Smith for- 
mulae 161, 



G{9 h 9 r , Vr) ^ H{6[ - n/2)H(6>. - 7 r/2)G , (9 i )G'(9 r ) 

where H is the Heaviside step function and G' is a function which is given 
by 

G(0) = - 7= 2 = = (10) 



1 + VI + 7 2 tan 2 9 

if the surface irregularities are distributed according to the Trowbridge- Reitz 



probability function 11 



3. Experimental results 

We observed the light reflected by samples of fluoropolymers for various 
incident and viewing angles, both in and out of the plane of incidence 0. 
Figure [2] shows the reflection distribution functions for three different sam- 
ples of polished PTFE and the copolymers PFA, FEP and ETFE, for the 
xenon scintillation light (175 nm), in the plane of incidence. The light was 
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Viewing Angle 6 r Viewing Angle 9 r Viewing Angle 9 r 



Figure 2: Reflectance distributions of various fluoropolymers samples for the xenon scintil- 
lation light, produced as indicated, plotted as a function of the viewing angle (in degrees), 
for various angles of incidence. Negative values of r mean that ip = tt in figure [IJ 

produced in a proportional counter filled with gaseous xenon. The experi- 
mental procedure and the chamber used for the measurements are described 
in detail in refs. [Ej and [3]. 

The reflectance distributions of figure [2] show clearly the presence of the 
three contributions: a diffuse lobe, a specular lobe and a specular spike. 
The data was fitted with the function g and the values of the four unknown 
parameters extracted, as shown in table [Q 

In figure [3] we compare the fits to a sample of unpolished PTFE using 
either one of the function distributions Ptr or Pqt (eqs. M and [9]). The 
agreement is clearly better if Ptr is used. Therefore, we used this function 
in all the fits shown in the figure [2j 

PTFE and PFA show a pronounced diffuse lobe, corresponding to an 
albedo between 0.6 to 0.7. The intensity of this component varies as expected 
with the angle of incidence, being suppressed at large angles in result of 
decreasing the intensity of the transmitted wave and consequently the sub- 
surface scattering. 
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None of the samples measured is a pure diffuser. The coherent spike 
increases with the angle as expected and it is especially visible at large angles 
of incidence. Its observed width is merely instrumental due to the aperture 
of the incident beam and of the aperture of the photo-detector used in the 
experiment. 

The effects of shadowing and masking reduce the specular lobe, specially 
at grazing angles, in which case the light is scattered towards elsewhere. The 
geometric factor G referred to above is relevant at those angles (for angles in 
excess of 80°) to reproduce the observed suppression. 

The total reflectance of the surface can be calculated for each angle of 
incidence by integrating the function q for all viewing directions: 

R{9i) = J J ^Q(9i,9 r ,Lp r )sm9 r d6 r dLp r (11) 

It is shown in figure H]for polished pressed PTFE as a function of the angle of 
incidence, The three reflection components that are shown remain fairly 
constant up to 6^ ~ 60°, the diffuse reflection dominating all along at 68%, 
whereas the specular lobe and specular spike amount to only 2.8% and 0.5% 
respectively. Above 70° the coherent spike is the main component of the 
specular reflection. 

This model can be also applied to calculate the reflectance of PTFE in 
contact with a liquid instead of vacuum/gas. Assuming that the albedo and 
refraction index of PTFE remain the same it is straightforward to conclude 
that in this case much more light will be reflected by the surface. The ex- 
pected reflectance predicted by this model for a PTFE liquid-xenon interface 
is plotted in figure [5j 

4. Conclusion 

The reflectance of polished fluoropolymers exhibit three main reflection 
components, a diffuse lobe, a specular lobe and a specular spike. The process 
can be modeled with only four free parameters and reproduces fairly well the 
details of the reflectance distribution measured in these materials. 

The data show that the specular lobe is best reproduced if the reflecting 
surface is approximated in the model by an ensemble of ellipsoid-shaped 
micro-surfaces. This conclusion comes at no surprise though, since these 
materials are granular by nature. 
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The relative intensity of the coherent specular spike appears to vary expo- 
nentially with the cosine of the angle of incidence, cos0j. This is in contrast 
to the prediction of the Beckmann-Spizzichino model of a dependency in 
cos 2 ^. This issue deserves additional studies in the future. 

In a PTFE vacuum interface the diffuse lobe is dominant for most of 
the incident angles. However, the total reflectance computed for PTFE in 
contact with liquid xenon shows a much more pronounced contribution from 
the specular reflection. The average reflectance is in agreement with previous 
estimates resulting in values of 90% for the total reflectance assumed to be 



independent of the incident angle, in this case [18j, [19 
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Figure 3: Comparison of the fits to non-polished PTFE using either the Trowbridge- Rcitz 
distribution function (red curves) or the Cook- Torrance distribution (blue curves). 




Angle of Incidence 9\ 



Figure 4: The reflectance of pressed polished PTFE as a function of the angle of incidence 
9\ for light of = 175 nm in vacuum. 
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Angle of Incidence 9\ 

Figure 5: The expected reflectance of pressed polished PTFE in contact with liquid xenon 
as a function of the angle of incidence 9i for light of = 175 11m (nx c = l-69[17j|). 
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